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We report a measurement of the fraction of events with a large pseudorapidity gap Ar; within the 
pseudorapidity region available to the proton dissociation products X inp + p^p + X. For a final 
state p of fractional momentum loss and 4-momentum transfer squared tp within 0.06 < < 0.09 
and \tp\ < 1.0 [0.2] GeV^ at ^ = 1800 [630] GeV, the fraction of events with Ar; > 3 is found to 
be 0.246 ± 0.001 (stat) ± 0.042 (syst) [0.184 ± 0.001 (stat) ± 0.043 (syst)]. Our resuhs are compared 
with gap fractions measured in minimum bias pp collisions and with theoretical expectations. 

PACS numbers: 13.85.Ni 



In a previous Letter [ij , we reported a measurement 
of the fraction of events with a central pseudorapidity 
gap At/ 2] produced in pp colhsions at -/s — 1800 
and 630 GeV. Here, we present results from a similar 
measurement performed in a sub-sample of pp events 
containing a leading (high longitudinal momentum) 
antiproton (Fig. 1). Large pseudorapidity gaps are 
presumed to be due to Pomeron {IP) exchange and are 
the signature for diffraction The process with a 
leading beam particle in the final state, which is kine- 
matically associated with an adjacent pseudorapidity 
gap, is known as single diffraction dissociation (SD), 
while that with a central gap as double diffraction dis- 
sociation (DD). The process in Fig. 1 is a combination 
of p-p SD and P-p DD and will be referred to in this 
paper as SDD. 

Low transverse momentum {pt) processes have tra- 
ditionally been treated theoretically in the framework 
of Regge theory l^. The introduction of a linear 
Pomeron {IP) trajectory, a{t) = a{Q) + a't, with inter- 
cept a{Q) = 1 -I- e > 1, enables the theory to correctly 
predict certain salient features of the high energy be- 
havior of hadronic interactions, such as the rise of the 
total cross section and the shrinking of the forward 
elastic scattering peak with increasing energy. How- 




FIG. 1: Schematic diagram and event topology in pseu- 
dorapidity space of a SDD (single diffraction plus gap) 
interaction, p + p p + GAPp + Afi + GAP + M2, with a 
leading outgoing antiproton of fractional momentum loss 
associated with a pseudorapidity gap Ar^p = In J^, and 

a gap within the region of 77 spanned by In s' — In s — In . 



ever, the success of the theory in describing diffrac- 
tion has been limited. While the shape of the SD and 
DD distributions as a function of A77 are correctly 
described, the normalization was found to be sup- 
pressed relative to the theoretical predictions by about 
an order of magnitude as the energy increases from 
~ 20 to 1800 GeV 0, S H. Proposals made to 
address this issue are divided into two general groups: 
(a) those based on "damping" of the SD cross sec- 
tion at small (anti)proton fractional momentum loss 
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^ or on a changing Pomeron intercept as a function 
of y/s or as a function of ^ 0, and (b) those in 
which only the overall normalization is suppressed as 
-/s increases [HSSEl- The models of group (a) can- 
not be applied to SDD. In a parton model approach 
developed by Bjorken to describe events with a 
large pseudorapidity gap between two jets produced 
by a colorless two-gluon exchange, a suppression of 
the overall normalization relative to the QCD calcu- 
lation was predicted due to additional partonic color 
exchanges in the same event which spoil the diffractive 
pseudorapidity gap signature. In this model, which 
belongs to group (b), the color exchanges would si- 
multaneously spoil all diffractive rapidity gaps in an 
event and therefore the ratio of two-gap to one-gap 
events in a given interaction would be unaffected. This 
approach should also hold for low pt diffractive pro- 
cesses. In this Letter, we examine this issue by study- 
ing p+p — *■ p-\-X with/without a gap A77 within the 77 
range of the system X in addition to the gap of nom- 
inal value Aryp = — In expected to be associated 
with the leading final state antiproton. 

Our study is based on our =1800 (630) GeV 
inclusive SD data described in Ref. [13 (El)- The 
events were collected in the 1995-96 Tevatron Run IC 
by triggering the Collider Detector at Fermilab (CDF) 
on an antiproton detected in a Roman Pot Spectrome- 
ter (RPS) [13. The average instantaneous luminosity 
during event collection was 0.2 x 10^° (1.5 x 10^°) cm~^ 
sec"^ at =1800 (630) GeV. The components of 
CDF 01 relevant to this study are the central track- 
ing chamber (CTC), the calorimeters, and two scintil- 
lation beam-beam counter (BBC) arrays. The CTC 
tracking efRciency varies from ^ 60% for pT — 300 
MeV to over 95% for pr > 400 MeV within \r]\ < 1.2, 
and falls monotonically beyond I77I = 1.2 approaching 
zero at I77I ^ 1.8. The calorimeters have projective 
tower geometry and cover the regions \ri\ < 1.1 (cen- 
tral), 1.1 < < 2.4 (plug), and 2.2 < I77I < 4.2 (for- 
ward). The Ar] x A(f) tower dimensions, where (j) is the 
azimuthal angle, are approximately 0.1 x 15° for the 
central and 0.1x5° for the plug and forward calorime- 
ters. The BBC arrays cover the region 3.2 < jryj < 5.9. 

The events were required to have a reconstructed 
RPS track of 0.06 < Ip < 0.09 and \tp\ < 1.0 [0.2] 
GeV^ at = 1800 [630] GeV, a hit on BBCp 
(proton-side BBC) to exclude DPE events, and no 
more than one reconstructed vertex within ±60 cm 
from the center of the detector along the beam di- 
rection. The vertex requirement was imposed to re- 
ject overlap events due to multiple interactions in the 
same beam-beam crossing, since additional interac- 
tions would most likely spoil the pseudorapidity gap 
signature of a diffractive event. 

At =1800 (630) GeV, the average value of 
0.075 of our data samples corresponds to IP-p collision 
energies of Vs' — \/^pS — 493 (173) GeV, at which 



the proton dissociation products cover the nominal 
?7-range from the maximum of 77 = In-^/s = 7.5 (6.5) 
down to 7/ = - In(^Vs) = -4.9 (-3.9). Thus, the CDF 
calorimeter coverage, \ri\ < 4.2, is well suited for the 
present study. 

Our analysis is similar to that used in evaluating the 
DD fraction in minimum bias events collected with a 
BBC coincidence trigger The method we use is 
based on the approximately flat dependence of the 
event rate on A77 expected for SDD events compared 
to the exponential dependence expected for the nor- 
mal SD events where rapidity gaps within the diffrac- 
tive cluster X are due to random multiplicity fluctu- 
ations. Thus, in a plot of event rate versus Ary, the 
SDD signal will appear as a flattenin g of an exponen- 
tially falling distribution at large A77 [ll| . In order to 
independently monitor detector effects in the positive 
and negative rj directions, we look for rjmax (Vmin), 
the 77 of the particle closest to 77 = in the pro- 
ton (antiproton) direction, and measure experimen- 
tal gaps overlapping 77 = 0, Ar]°^p = r]max - Vmin 
(see Fig. 1). For this purpose, a particle is defined 
as a reconstructed track in the CTC, a calorimeter 
tower with energy above a given threshold, or a BBC 
hit. The tower energy thresholds used, chosen to lie 
comfortably above noise level, are Et = 0.2 GeV for 
the central and plug and E = 1 GeV for the forward 
calorimeters. The calorimeter noise was measured us- 
ing beam-beam crossing events with no reconstructed 
vertex. At the calorimeter interfaces near |7/| ~ 0, 
1.1 and ~ 2.4, where the noise level was found to be 
higher, we use higher thresholds of up to 0.3 GeV. 
The average number of calorimeter towers per unit 
A77 with EJ^f'^'^^ above threshold in an event is ~ 0.07, 
small compared to the corresponding average particle 
density of ~ 3 in the data. The fraction of SDD to 
total number of events based on Arj^^p is obtained di- 
rectly from the data and corrected for (a) contamina- 
tion from SD events, (b) the effect of the unobserved 
(below threshold) particles, and (c) the triggering effi- 
ciency (acceptance) of BBCp for SD and SDD events. 
These corrections are made using a hadron-level MC 
simulation. 

The MC generation of SD events is described in 
Ref. Q. For SDD events, a diffractive p and a cluster 
of = are generated as for a SD interaction, and 
a DD interaction is assumed to takeplace in the P-p 
colhsion, which is treated as in Ref. jlj and boosted to 
the lab frame. The same thresholds are used for parti- 
cles in the MC as for towers in the data, after dividing 
the generated particle Et by an 77-dependent calibra- 
tion coefficient of average value ~ 1.6 representing 
the ratio of true to measured calorimeter energy. The 
MC generator includes the calorimeter noise, and for 
charged particles it is followed by a detector simula- 
tion. 

Figure 2 shows lego histograms of events versus 
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Vmax and —rjjnin for (a) data and (b) MC generated 
events, as well as MC events for (c) only SD and (d) 
only SDD at ^/s = 1800 GeV. Similar results are ob- 
tained at ^/s = 630 GeV. The observed structure in 
the distributions along rj^y^ciximm) is caused by the vari- 
ation of the tower energy thresholds with \r]\. The bins 
at Irjmaximin) I = 3.3 Contain all events within the BBC 

range of 3.2 < \T]max(rmn) \ < 5.9. 

Vs=1800 GeV 




FIG. 2: The number of events as a function of rjmax and 
—rjmin, the r] of the track or hit tower closest to = 
in the proton and antiproton direction, respectively, at 
^/s = 1800 GeV: (a) data; (b) Monte Carlo simulation; 
(c,d) the individual contributions of Monte Carlo gener- 
ated SD and SDD (SD -I- gap) events. The MC distribu- 
tions are normalized by a two-component fit to the data 
using distributions (c) and (d). 

Figure 3 presents the number of events as a function 
of Ar^g^p for the 1800 GeV data (points) and for a fit 
to the data using a mixture of MC generated SD and 
SDD contributions (solid histogram). The SD contri- 
bution (dashed histogram) exhibits an approximately 
exponential fall with increasing Ar^^^p, as expected. 
The region of Arj'^^p > 3 is dominated by the SDD 
signal and is used to extract the gap fraction (ratio of 
SDD to total number of events). The approximately 
flat behaviour expected for the SDD distribution in 
this region is modulated by the ry-dependent tower 
thresholds used, causing the observed bumps and dips, 
and by the BBCp acceptance for SDD events, which 
decreases with increasing A77°^p. The MC simulation 
reproduces these features of the data. 

At = 1800 [630] GeV, the fraction of events with 
AryO^p > 3 is (15.9 ± 0.1)% [(17.5 ± 0.2)%], of which 
the contribution of background SD events, estimated 
using the MC simulation, is 1.2% [2.4%]. The quoted 
errors are statistical. The amount of SD background 
in the region Ar]'^ > 3 depends on the tower energy 
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max imin 

FIG. 3: The number of events as a function of Arj^^p = 
rimax — rjmin for data at -/s = 1800 GeV (points), for SDD 
(SD 4- gap) MC generated events (solid line), and for only 
SD MC events (dashed line). 

calibration coefficients and thereby on the calorimeter 
tower energy thresholds used in the MC. For example, 
increasing these thresholds has the effect of decreasing 
the multiplicity in the MC generated events, result- 
ing in larger rapidity gaps and hence larger SD back- 
ground in the region of Arf^^^ > 3. The systematic 
uncertainty in the background is estimated by raising 
(lowering) the tower thresholds in the MC by a fac- 
tor of 1.25 (evaluated from a multiplicity uncertainty 
of ±10%), which increases (decreases) the background 
by a factor of 1.6. 

The correction factors needed to account for the 
effect of unobserved particles and thus convert the 
measured gap fractions to gap fractions correspond- 
ing to the nominal gap definition '3| , Arp = In j^j-i'^^-i 

(In . ; < 0, i — 1,2), where sq = 1 GeV^, were 

V S So 

evaluated using the SDD Monte Carlo simulation and 
found to be 0.81 [0.73] at = 1800 [630] GeV. 
Systematic errors are obtained by varying the tower 
energy thresholds by ±25%. These errors are corre- 
lated with the errors in the SD background contamina- 
tion and therefore a combined systematic uncertaity 
of 15% [23%] was evaluated for both these effects and 
applied to the extracted nominal gap fractions. 

The BBCp acceptance, evaluated from the SDD 
MC simulation, is 0.68±0.06(syst) [0.81 ±0.04(syst)], 
where the error is due to a 20% systematic uncertainty 
assigned to the fraction of the diffraction dissociation 
mass clusters which do not trigger the BBCp. For 
SD events, the BBC acceptance is 0.98 ± O.Ol(syst) 
[0.98 ± O.Ol(syst)]. Including all systematic errors, 
the acceptance-corrected SDD fractions for nominal 
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gaps At]° > 3 are 0.174 ± 0.001 (stat) ± 0.030 (syst) 
[0. 138 ± 0.001 (stat) ± 0.032 (syst)] at = 1800 [630] 
GeV. 

The Ajf > 3 fractions are extrapolated to all SDD 
gaps of A77 > 3 using the shape of the gap distribu- 
tion of Eq. 1, which is based on Regge theory and 
factorization. This equation, which was used in the 
MC simulation, is obtained from the equation for SD 
by replacing the IP-p total cross section factor with 
the P-p DD factor k{...} (see Ref. Q). 



dtpdtd^pdArjdric 



m 



n 2 



[Q(t)-l]Ar, 



/3^(0)(- 



(1) 



Here, rjc is the center of the gap At], /3(0) the IP-p 
coupling, K the ratio of the triple-Pomeron to the IP- 



p couplings, and In = In ^ — In ^i- — Ar/ the rapidity 

space occupied by particles; \/s" will be referred to be- 
low as "diffractive sub-energy" . For numerical evalu- 
ations we use Q e = 0.104 ±0.002, a' = 0.25 GeV-^, 
K = 0.17, P{tp) = 6.57 GeV"^ [4.1 mb*] x Fi{tp), 
where Fi(tp) is the nucleon form factor. The variable 
t is not measured and therefore is integrated over in 
the calculations. Owing to the increased phase space 
resulting from releasing the requirement that the ra- 
pidity gap overlap = 0, the Arj > 3 fractions are 
found to be larger than the Ayy^ > 3 ones by a factor 
of 1.44 [1.40] at Vs = 1800 [630] GeV. The evaluation 
of this factor is performed analytically and therefore 
no error is assigned to it; the effect of the uncertainty 
in the parameter e, which controls the shape of the 
Ar] distribution in Eq. 1, is < 1%. 

Our results for the ratio Rg^,^ of the number of 
events with a gap of A-q > 3 to the total number 
of SD events are 0.246 ± 0.001 (stat) ± 0.042 (syst) 
[0.184 ±0.001 (stat) ±0.043 (syst)] at ^/i = 1800 [630] 
GeV. These ratios are plotted in Fig. 4 at Vs' — 493 
[173] GeV, the average value of the diffractive mass 
Mx, and compared with double diffractive to total 
cross section ratios — /ar, where is ob- 
tained from [3] and ctt is set to the Pomeron exchange 

contribution, /3^(0) (yj^ ■, to conform with the defini- 
tion of R%%^ ■ The vertical error bars are mainly due 
to systematic effects, which are correlated among all 
points. The dashed lines represent predictions based 
on Regge theory and factorization normalized to the 
SD cross section at y/s =22 GeV (see ^3)- The solid 
lines are predictions from the "renormalized gap prob- 
ability" model , in which the Regge cross section is 
factorized into two parts, one representing the pp total 
cross section at the diffractive sub-energy multiplied 
by k", where n is the number of gaps, and a factor 
interpreted as the gap probability distribution normal- 
ized to unity over all available phase space. The bands 
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FIG. 4: Ratios of SDD (single-difltraction plus gap) to 
single-difFractive rates (filled circles) and double-diffractive 
to total cross sections (open circles) as a function of the 
collision energy of the sub-process, Pomeron-proton and 
pp, respectively. The uncertainties are mainly due to sys- 
tematic effects, which are highly correlated among all four 
data points. The dashed lines are predictions from Regge 
theory and the solid lines from the renormalized gap prob- 
ability model [lo|. 



around the solid lines represent a 10% uncertainty due 
to the factor k The data are in good agreement 
with the renormalized gap model predictions. 

In summary, we have measured the fraction of 
events with a pseudorapidity gap Ar] > 3 within 
the diffractive cluster X of the process pp pX 
and found it to be 0.246 ± 0.001 (stat) ± 0.042 (syst) 
[0.184 ± 0.001 (stat) ± 0.043 (syst)] for 0.06 < ^,5 < 
0.09 and \tp\ 



< 1.0 [0.2] GeV^ at ^ =1800 [630] 
GeV. These values are higher than expectations from 
double-diffractive fractions in minimum bias events, 
lower than expectations from Regge theory and fac- 
torization, and in good agreement with predictions 
based on the renormalized gap probability model 10]. 
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